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Abstract 

Current  density  and  potential  distribution  measurements  were  conducted  using  a  segmented  current  collector  and  flow  field  setup  on 
membrane  electrode  assemblies  prepared  with  segmented  and  unsegmented  electrodes  made  from  two  different  types  of  commercially 
available  gas  diffusion  layers.  Both  galvanostatic  and  potentio static  discharge  modes  were  employed.  Irrespective  of  the  type  of  gas  diffusion 
layer,  when  a  common  electrode  was  employed,  significant  performance  variations  were  encountered  between  current  collector  segments 
in  the  constant  voltage  mode,  while  the  segment  to  segment  variations  were  minimal  in  the  constant  current  mode.  Both  types  of  discharge 
modes  showed  negligible  variations  between  segments  in  the  case  of  segmented  electrode.  A  simple  mathematical  model  was  developed 
to  assist  in  the  interpretation  of  the  experimental  results.  The  differences  in  contact  resistances  between  the  current  collectors  and  the  gas 
diffusion  layer,  especially  on  the  cathode  side  have  been  identified  as  the  primary  reason  for  the  experimentally  observed  behavior.  Based 
on  the  results  presented  here,  segmenting  the  electrode  along  with  the  current  collector  is  recommended  for  current  distribution  studies. 
When  using  a  common  electrode,  only  the  galvanostatic  mode  is  preferred  to  minimize  contact  artifacts. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cells  are  fast  becoming  viable  alternative  energy 
conversion  devices  of  the  21st  century.  Their  high  effi¬ 
ciency,  simplicity  in  design  and  operation  and  pollution  free 
characteristics  make  them  an  attractive  option  for  terrestrial 
applications.  Of  the  various  types  of  fuel  cells  available, 
proton  exchange  membrane  (PEM)  fuel  cells  are  consid¬ 
ered  to  be  most  suitable  for  transportation  and  portable 
applications  due  to  attractive  features  like  low  operation 
temperature,  high  energy  density  and  efficiency.  However, 
their  performance  needs  to  be  further  optimized  [1]  to  be 
cost  competitive  with  current  energy  conversion  devices 
like  the  internal  combustion  engine  or  batteries.  Significant 
strides  have  been  made  towards  addressing  this  goal  of 
performance  improvement  and  cost  optimization  through 
development  of  better  membranes,  improved  catalyst  layer 
fabrication  techniques  leading  to  better  catalyst  utilization 
[1]  and  better  flow  field  designs  that  enhance  reactant  and 
product  transport  within  the  fuel  cell  [1,2]. 
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Performance  improvement  in  a  typical  PEM  fuel  cell  di¬ 
rectly  translates  to  achieving  higher  average  current  per  unit 
area  (otherwise  known  as  current  density)  at  any  given  oper¬ 
ating  cell  voltage  or  vice  versa.  On  the  other  hand,  uniformity 
of  current  density  distribution  across  the  entire  active  area 
is  crucial  for  performance  optimization.  The  local  current 
density  distribution  within  a  PEM  fuel  cell  is  a  function  of 
various  factors  like  local  membrane  hydration  state,  reactant 
and  product  concentration,  temperature,  etc.  Over  the  past 
decade,  experimental  complications  associated  with  PEM 
fuel  cells  due  to  their  inherent  small  geometries  and  extreme 
aspect  ratios,  prompted  researchers  to  develop  representative 
first  principles-based  mathematical  models  [3-18],  to  gain 
qualitative  information  on  membrane  hydration,  temperature 
and  species  concentration  distribution.  Pioneering  work  in 
fuel  cell  modeling  were  usually  one  dimensional,  represent¬ 
ing  the  direction  normal  to  the  reactive  catalyst  surface,  and 
accounted  only  for  gaseous  phase  to  avoid  the  complexi¬ 
ties  involved  in  multi-dimensional  modeling  of  multi-phase 
flow  in  porous  media  [3-6].  Though  these  models  provide 
excellent  qualitative  information,  to  achieve  representation 
of  ‘real  life’  situations,  researchers  have  been  looking  into 
multi-dimensional,  multi-phase  models  [7-18]  to  get  a  more 
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Nomenclature 

equilibrium  potential  (V) 

k 

apparent  exchange  current  density 
(A/cm2) 

^APP 

applied  current  (A/cm2) 

^CON,A 

areal  contact  resistance  between  current 
collector  A  and  gas  diffusion  layer 
(Q  cm2) 

^CON,B 

areal  contact  resistance  between  current 
collector  B  and  GDL  (£2  cm2) 

Lapp 

applied  voltage  (V) 

AV 

total  voltage  drop  from  the  sensing  point 
to  the  reactive  interface  (V) 

AVcon 

voltage  drop  due  to  contact  resistance  (V) 

AVgdl 

voltage  drop  across  the  gas  diffusion 
layer  (V) 

AVk 

voltage  drop  due  to  activation  (V) 

Greek  letters 

0A 

empirical  anode  kinetic  slope  (V/A/cm2) 

OC 

cathode  Tafel  slope  (V/decade) 

^CC 

bulk  conductivity  of  POCO  graphite 
current  collector  (S/cm) 

G  X,GDL 

through-plane  conductivity  of  gas 
diffusion  layer  (S/cm) 

&Y,  GDL 

in-plane  conductivity  of  gas  diffusion 
layer  (S/cm) 

0 

potential  (V) 

accurate  picture  of  the  various  transport  and  kinetic  phenom¬ 
ena  in  PEM  fuel  cells  employing  different  flow  distribution 
strategies.  While  these  models  provide  excellent  qualitative 
information,  experimental  data  on  the  local  current  density 
distribution  in  operating  PEM  fuel  cells  is  essential  to  vali¬ 
date  and  verify  the  model  predictions  and  to  accurately  esti¬ 
mate  the  various  kinetic  and  transport  parameters  to  develop 
these  models  into  practical  design  tools. 


2.  Background 

In  this  section  available  literature  on  attempts  to  obtain 
such  local  current  density  distribution  data  is  reviewed 
briefly.  Cleghorn  et  al.  [19]  conducted  some  pioneering  cur¬ 
rent  density  distribution  measurements  on  typical  lab  scale 
PEM  fuel  cell  setups.  The  authors  used  printed  circuit  board 
(PCB)  technology  to  create  a  segmented  current  collector 
and  flow  field  that  was  used  on  the  anode  with  segmented 
gas  diffusion  layer  and  catalyst  layer,  while  the  cathode 
employed  a  regular  unsegmented  electrode,  (i.e.  diffusion 
and  catalyst  layer)  current  collector  and  flow  field.  The 
effects  of  anode  and  cathode  stream  humidity  and  stoichio¬ 
metric  flow  rate  of  air  on  the  steady-state  current  density 


distribution  were  studied.  The  authors  employed  a  com¬ 
bination  of  two  load  units  and  a  specially  designed  patch 
board  that  acted  as  a  multiplexer  to  control  the  voltage  at 
the  various  segments.  Stumper  et  al.  [20]  analyzed  three 
methods  for  current  density  distribution  mapping  namely 
the  partial  membrane  electrode  assembly  (MEA)  approach, 
the  subcell  technique  and  the  current  mapping  technique. 
The  first  approach  involves  the  use  of  several  different 
MEAs  with  a  catalyzed  active  area  of  varying  fractions  of 
the  total  flow  field  area.  In  the  subcell  approach  the  authors 
used  a  number  of  ‘  subcells  ’  at  various  locations  along  the 
gas  flow  channel  that  were  electrically  insulated  from  the 
main  active  MEA  and  controlled  by  a  separate  load.  In  the 
third  approach,  a  network  of  passive  graphite  resistors  were 
placed  between  the  flow  field  plate  and  the  current  collect¬ 
ing  bus  plate,  while  the  potential  drop  across  these  resistors 
were  monitored  to  establish  the  current  flowing  through 
them.  Wieser  et  al.  [21]  proposed  the  use  of  an  array  of 
Hall  sensors  with  a  segmented  current  collector  and  flow 
field  for  local  current  density  measurements.  The  authors 
used  unsegmented  electrodes  in  their  study. 

While  these  techniques  are  no  doubt  quite  innovative,  they 
are  not  without  disadvantages.  The  partial  MEA  approach 
does  not  provide  sufficient  spatial  resolution  and  significant 
errors  can  arise  due  to  inherent  variations  in  electrical,  trans¬ 
port  and  kinetic  properties  between  different  MEAs.  The 
‘subcelT  approach  is  plagued  by  the  difficulty  in  properly 
isolating  the  4  subcells  ’  from  the  main  electrode  and  achiev¬ 
ing  perfect  alignment  of  the  anode  and  cathode  sides.  The  use 
of  Hall  sensors  can  significantly  complicate  the  experimental 
setup  making  experimentation  expensive  and  tedious.  More¬ 
over,  interference  from  neighboring  segments  is  also  possi¬ 
ble  with  the  use  of  Hall  sensors.  A  major  cause  of  concern 
with  the  current  mapping  technique  using  passive  resistors 
is  the  lateral  in-plane  currents  through  the  flow  field  plate 
that  could  lead  to  very  low  spatial  resolution  due  to  current 
spreading.  The  issue  of  significant  lateral  currents  can  be 
further  complicated  by  differences  in  the  contact  resistance 
between  the  various  resistors  and  the  flow  field  plate.  When 
using  a  single  potentiostat,  if  the  voltage  sensing  and  con¬ 
trol  point  is  after  the  resistors  the  potential  at  the  electrode 
will  not  be  uniform  because  of  the  differences  in  the  current 
flowing  through  the  resistors.  Though  Cleghorn  et  al.’s  [19] 
approach  avoids  these  complications,  the  use  of  a  couple  of 
load  units  with  a  multiplexer  (i.e.  voltage  control  is  quickly 
switched  from  one  electrode  to  another  by  the  multiplexer) 
only  allows  analysis  of  steady-state  behavior.  Switching  be¬ 
tween  electrodes  gives  rise  to  temporal  double  layer  charg¬ 
ing  and  discharging  currents  and  could  induce  undesired 
transient  artifacts.  Barring  the  partial  MEA  approach,  the 
problem  associated  with  the  use  of  a  single  potentiostat  is 
common  to  all  the  other  methods  mentioned  above. 

Natarajan  and  Nguyen  [22]  presented  some  qualitative  re¬ 
sults  on  hydrogen  starvation  effects  using  segmented  current 
collector  and  segmented  electrode  (diffusion  and  catalyst 
layers)  setup  on  the  anode  side  while  using  a  single  common 
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electrode  and  current  collector  on  the  cathode  side.  The  au¬ 
thors  employed  six  segments  along  a  single  straight  channel 
(corresponding  to  the  conventional  or  serpentine  flow  field 
design)  that  were  individually  controlled  by  a  multi-channel 
potentiostat/galvanostat  to  avoid  the  possible  pit  falls  associ¬ 
ated  with  the  use  of  a  single  load  and  a  multiplexer.  Further 
details  regarding  the  experimental  setup  are  discussed  in  the 
following  sections.  Brett  et  al.  [23]  adopted  the  PCB  tech¬ 
nique  introduced  by  Cleghorn  et  al.  [19]  to  a  single  straight 
channel  cathode  where  each  of  the  current  collector  ribs  (10 
in  all)  on  the  PCB  board  was  controlled  by  an  individual 
load.  They  used  a  typical  multi-channel  anode  with  hydro¬ 
gen  flowing  perpendicular  (cross-flow)  to  the  air  stream.  The 
electrodes  on  both  sides  of  the  MEA  were  not  segmented. 
The  authors  have  reported  both  steady-state  and  transient 
results  in  their  publication. 

Recently  Mench  and  Wang  [24,25]  published  some  inter¬ 
esting  steady- state  and  transient  current  density  distribution 
data  in  direct  methanol  fuel  cells  [24]  and  PEM  fuel  cells 
[25].  They  used  specially  built  segmented  flow  fields  with 
serpentine  flow  channels  for  both  the  anode  and  cathode 
sides,  where  the  individual  segments  were  separately  con¬ 
trolled  by  a  multi-channel  potentiostat/galvanostat.  The  seg¬ 
mented  flow  fields  were  fabricated  by  embedding  gold  plated 
stainless  steel  ribs  in  a  polycarbonate  block.  The  authors  em¬ 
ployed  commercially  available  MEAs  where  the  electrodes 
were  not  segmented.  Nopenen  et  al.  [26]  and  Hottinen  et  al. 
[27]  adopted  a  similar  approach  to  study  steady- state  and 
transient  current  density  distributions  in  free-breathing  PEM 
fuel  cells. 

Modifications  like  the  use  of  segmented  current  collec¬ 
tors  and  individual  control  of  segments  using  multiple  po- 
tentiostats  to  the  innovative  techniques  proposed  by  Stumper 
et  al.  [20]  and  Cleghorn  et  al.  [19]  eliminate  problems  like 
in-plane  current  spreading  within  the  flow  field  plate  or  un¬ 
desired  transients  caused  by  the  multiplexer.  However,  apart 
from  these  advances  one  other  major  question  to  be  ad¬ 
dressed  in  current  density  distribution  studies  is  whether  the 
gas  diffusion  layer  (GDL)  and  the  catalyst  layer  (i.e.  the 
electrode)  need  to  be  segmented  along  with  the  current  col¬ 
lector. 

Current  density  distribution  studies  are  aimed  at  obtain¬ 
ing  information  on  the  local  reaction  rate  distribution  within 
the  catalyst  layer  which  is  a  strong  function  of  the  local 
electronic  and  ionic  potentials,  reactant  and  product  con¬ 
centration  and  membrane  hydration  state.  Performance  im¬ 
provement  and  optimization  of  fuel  cells  directly  relates  to 
improving  the  electrochemical  reaction  rate  within  the  cat¬ 
alyst  layer.  Hence  an  accurate  picture  of  the  local  current 
density  distribution  as  it  emanates  from  the  catalyst  layer 
is  of  paramount  interest  to  researchers  from  a  fundamental 
point  of  view.  In  the  case  of  a  common  electrode  the  ques¬ 
tion  of  whether  one  sees  the  same  current  distribution,  as  in 
the  catalyst  layer,  at  the  current  collectors  depends  on  (1) 
the  ratio  of  the  normal  distance  from  the  current  collector  to 
the  catalyst  layer  and  electrode  width  along  the  channel,  (2) 


the  ratio  of  in-plane  and  through-plane  area,  (3)  the  ratio  of 
in-plane  and  through-plane  conductivities  of  the  GDL  mate¬ 
rial  and  finally,  (4)  the  relative  contact  resistances  between 
the  GDL  and  the  current  collectors.  The  above  mentioned 
parameters  determine  the  extent  of  current  spreading  within 
the  GDL.  Mench  and  Wang  [24,25]  implemented  elaborate 
experimental  techniques  like  the  use  of  pressure  indicating 
film  to  minimize  contact  resistance  and  assumed  minimal 
current  spreading  or  ‘crosstalk’  in  their  experiments  with 
common  GDLs.  Another  concern  with  the  use  of  a  common 
GDL  is  determining  the  active  area  associated  with  each 
current  collector.  Nopenen  et  al.  [26]  tried  to  address  this 
by  developing  a  simple  model  to  estimate  the  errors  arising 
from  assigning  equal  areas  to  all  the  current  collectors. 

Experimental  verification  and  quantification  of  the  exis¬ 
tence  and  extent  of  current  dispersion  within  a  common  GDL 
in  contact  with  multiple  isolated  current  collectors  are  essen¬ 
tial  and  is  the  focus  of  this  work.  This  paper  also  emphasizes 
the  need  for  segmenting  the  electrode  in  order  to  minimize 
the  effect  of  variations  in  contact  resistances  between  the 
electrode  and  the  current  collectors.  A  simple  mathematical 
model  is  also  presented  here  to  assist  in  the  interpretation  of 
the  experimental  results.  The  need  for  segmenting  the  elec¬ 
trode  to  correctly  measure  the  true  current  density  distribu¬ 
tion  within  the  catalyst  layer  is  also  emphasized.  Finally,  the 
significance  of  these  results  in  terms  of  realistic  PEM  fuel 
cell  model  development  is  also  discussed. 

3.  Experimental 

3.1.  Segmented  current  collector  and  flow  field 

Fig.  1  is  a  schematic  of  the  top  and  side  view  of  the 
current  collector  and  flow  field  block  used  in  this  work. 
Using  a  computer  numerically  controlled  (CNC)  milling 
machine  six  rectangular  grooves  with  precise  dimensions  of 
50  mm  (length)  x  10  mm  (width)  x  4.85  mm  (depth)  were 
machined  into  a  140  mm  x  70  mm  block  of  acrylic.  Bolt 
holes  and  pipe  threaded  gas  holes  were  also  machined  into 
the  block.  Strips  of  the  above  mentioned  dimensions  were 
also  machined  out  of  a  5  mm  thick  POCO  graphite  plate 
using  a  CNC  milling  machine.  The  graphite  plates  were 
then  gently  pressed  into  these  slots  using  a  precision  vice 
and  quick  setting  ‘super  glue’  was  wicked  into  the  gaps 
between  the  grooves  and  graphite.  The  protruding  graphite 
surfaces  were  then  machined  flat  using  a  sharp  carbide  ‘fly 
cutter’.  After  polishing  the  surface  on  a  fine  sand-paper,  a 
single  gas  channel  was  machined  connecting  the  inlet  and 
exit  holes  as  shown  in  the  figure.  The  acrylic  block  was 
2.5  cm  thick  for  the  anode  while  the  cathode  graphite  strips 
were  placed  in  a  1.25  cm  thick  plastic  block.  A  blank  2.5  cm 
thick  acrylic  block  and  an  aluminum  heating  block  were 
also  machined  to  contain  pipe  threaded  gas  holes  and  bolt 
slots.  During  cell  assembly  the  aluminum  heating  block 
was  placed  on  the  cathode  side  with  O-rings  on  the  gas 
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Bolt  holes 


2  mm 


Fig.  1.  Schematic  of  the  current  collector  and  flow  field  block:  (a)  top  view;  (b)  side  view. 


holes  between  the  blank  plate  and  the  thinner  acrylic  base 
plate  to  achieve  better  heat  transfer.  The  major  planar  faces 
of  all  the  above  mentioned  parts  of  the  fuel  cell  assembly 
were  machine  finished  and  polished  and  were  found  to  be 
even  within  12.5  \xm  using  a  high  precision  dial  indicator. 

3.2.  ME  A  fabrication 

Electrodes  fabricated  by  catalyst  coating  on  two  different 
types  of  gas  diffusion  layers,  namely  SIGRACET®  GDL 
30  BC  from  SGL-CARBON  Inc.  and  Toray®  carbon  pa¬ 
per  TGPH-120  from  ETEK  Inc.  were  obtained  from  TVN 
Systems  Inc.  The  carbon  papers  were  330  and  350  jam 
thick,  respectively.  The  SIGRACET®  diffusion  layer  had  a 
Teflon®  content  of  5  wt.%  while  the  Toray®  contained  no 
added  Teflon®.  The  catalyst  loading  on  both  types  of  GDL 
were  about  0.35gPt/cm2.  A  schematic  of  the  MEAs  with 
segmented  or  unsegmented  electrodes  on  both  sides  are 
provided  in  Fig.  2.  Electrode  pieces  (7  mm  x  10  mm)  were 
cut  from  the  procured  stock  along  with  gaskets  with  slots 
for  the  active  electrode  and  holes  for  bolts,  using  precise 
paper  masks  with  dimensions  marked  on  it  using  AUTO¬ 
CAD.  Using  the  gaskets  as  a  frame  to  accurately  arrange 
the  electrodes,  the  electrodes  were  hot  pressed  to  Nation® 
112  membranes  at  135  °C  and  65psi  (based  on  total  area 
including  gaskets).  A  combination  of  Teflon®  and  silicone 
gaskets  were  used  on  both  sides  of  the  MEA.  Though  the 


fuel  cell  assembly  was  capable  of  handling  six  segments, 
MEAs  were  prepared  to  correspond  to  four  current  collector 
segments  for  the  sake  of  clarity  of  experimental  data.  More¬ 
over,  it  should  be  noted  that  the  common  MEA  had  a  slightly 
greater  active  area  (ca.  15%)  than  the  segmented  MEA  due 
to  the  region  between  the  current  collectors.  It  should  be 
noted  that  in  the  segmented  MEA,  the  electrodes  on  anode 
and  cathode  were  segmented,  while  a  common  electrode 
was  used  on  both  sides  in  case  of  the  common  MEA. 

3.3.  Test  procedure 

The  machined  parts  of  the  fuel  cell  and  the  fabricated 
MEAs  were  carefully  assembled  to  ensure  proper  alignment 
between  the  current  collector  and  electrode  edges.  The  cell 
was  assembled  in  a  press  at  the  same  pressure  as  the  hot 
pressing  step  of  65  psi  based  on  the  entire  cell  area  to  avoid 
any  distortion  to  the  MEA  due  to  excess  compression.  The 
cell  was  tested  for  gas  leak  and  cross-over  before  actual 
testing.  Hydrogen  was  provided  to  the  anode  (thick  plate) 
through  a  bottle  of  de-ionized  water  held  at  60  or  70  °C 
while  oxygen  was  supplied  to  the  cathode  sparged  at  room 
temperature.  Hydrogen  and  oxygen  flow  rates  employed  in 
all  the  experiments  were  about  51.5  and  25.8  cm3 /min,  re¬ 
spectively.  These  numbers  translates  to  2  A/cm2  for  common 
MEA  with  an  active  area  of  3.22  cm2,  while  for  the  seg¬ 
mented  MEA  with  an  active  area  of  2.8  cm2  they  correspond 
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(a)  \ 

Teflon  /  silicone  gasket 

/ 


Fig.  2.  Schematic  of  the  membrane  electrode  assembly:  (a)  segmented; 
(b)  unsegmented. 

to  2.3  A/cm2.  The  four  current  collector  segments  that  were 
in  contact  with  the  electrode  were  each  controlled  indepen¬ 
dently  using  a  multi-channel  potentiostat/galvanostat  from 
Arbin  Systems  through  current  and  voltage  leads  connected 
directly  to  the  graphite  segments  protruding  on  one  side 
(see  Fig.  1).  The  ME  As  were  subjected  to  multiple  cycles 
alternating  between  constant  current  and  voltage  staircases 
at  30  °C.  The  staircase  limits  for  the  potentiostatic  staircases 
were  chosen  to  be  open  circuit  voltage  and  0.7  V  while  the 
galvanostatic  staircases  were  cycled  between  open  circuit 
and  100mA  per  segment  for  the  SIGRACET®  diffusion 
layer  and  50  mA  per  segment  for  the  Toray®  GDL.  These 
values  for  the  cycle  limits  were  so  chosen  to  limit  mass 
transfer  effects  caused  by  liquid  water  accumulation  on 
the  current  density  distribution.  During  these  cycles  the 
segments  were  held  at  each  current  or  voltage  for  3  min.  Fi¬ 
nally,  once  the  MEAs  were  sufficiently  massaged  (evaluated 
based  on  reproducibility  between  cycles)  the  data  from  the 
last  galvanostatic  and  potentiostatic  cycles  were  chosen  for 
analysis. 

4.  Results  and  discussion 

Fig.  3  provides  the  experimental  data  in  the  form  of  po¬ 
larization  curves  based  on  the  current  and  voltage  responses 
from  the  four  segmented  current  collectors,  using  a  ME  A 


Abs.  value  of  segment  current  density  (A/cm2) 


Fig.  3.  Galvanostatic  and  potentiostatic  segment  polarization  curves  using 
a  segmented  MEA  with  SIGRACET®  GDL.  H2  flow  rate  =  2.3  A/cm2,  O2 
flow  rate  =  2.3  A/cm2,  H2  humidifier  temperature  =  70  °C,  O2  humidifier 
temperature  =  25  °C  and  cell  temperature  =  30  °C. 

with  segmented  electrodes  fabricated  from  a  SIGRACET® 
gas  diffusion  layer.  The  solid  lines  represent  data  obtained 
from  galvanostatic  discharge,  while  the  dashed  line  indicates 
potentiostatic  discharge.  Similar  results  for  a  MEA  with  a 
common  (unsegmented)  electrode  using  the  same  experi¬ 
mental  setup  and  operating  conditions  are  provided  in  Fig.  4. 
The  current  densities  for  the  common  MEA  were  calculated 
assuming  that  the  area  above  the  plastic  separator  is  evenly 
divided  between  adjacent  current  collectors.  Fig.  5  provides 
the  total  current  from  the  entire  cell  (sum  of  the  currents 
from  each  collector  segment)  verses  the  average  voltage  at 
the  current  collectors.  Experimental  data  obtained  from  the 
MEAs  prepared  from  Toray®  diffusion  layers  are  provided 
in  the  latter  sections. 

Fig.  3  shows  that  in  the  case  of  the  segmented  MEA,  for 
any  given  segment  the  polarization  curves  obtained  by  both 
types  of  discharges  are  very  close  to  each  other  and  that  the 
segment  to  segment  performances  are  also  uniform  (a  spread 


Abs.  value  of  segment  current  density  (A/cm2) 


Fig.  4.  Galvanostatic  and  potentiostatic  segment  polarization  curves  using 
a  common  MEA  with  SIGRACET®  GDL.  H2  flow  rate  =  2  A/cm2,  O2 
flow  rate  =  2A/cm2,  H2  humidifier  temperature  =  70  °C,  O2  humidifier 
temperature  =  25  °C  and  cell  temperature  =  30  °C. 
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Abs.  value  of  total  current  (A) 

Fig.  5.  Comparison  of  total  cell  current  of  common  and  segmented  ME  As 
under  galvanostatic  and  potentiostatic  discharge  with  SIGRACET®  GDL. 
H2  flow  rate  =  2  or  2.3  A/cm2,  O2  flow  rate  =  2  or  2.3  A/cm2,  H2 
humidifier  temperature  =  70  °C,  O2  humidifier  temperature  =  25  °C  and 
cell  temperature  =  30  °C. 

within  5-7  mA  at  any  given  operating  voltage  or  vice  versa). 
On  the  other  hand,  Fig.  4  indicates  that  in  a  common  (unseg¬ 
mented)  MEA,  while  the  polarization  curves  from  the  four 
collectors  are  quite  close  to  each  other  in  the  constant  cur¬ 
rent  discharge  mode  (within  5  mV  at  100mA),  there  exists 
a  significant  spread  in  performance  from  segment  to  seg¬ 
ment  under  constant  voltage  mode  (40  mA  at  0.725  V).  The 
difference  between  the  segments  under  potentiostatic  mode 
seems  to  increase  with  decreasing  cell  voltage.  It  should  be 
noted  here  that  the  gas  flow  rates  in  these  experiments  were 
chosen  such  that  there  were  negligible  concentration  vari¬ 
ations  along  the  channel.  Moreover,  on  comparing  Figs.  3 
and  4  it  is  evident  that  the  polarization  curves  obtained  in 
the  galvanostatic  or  constant  current  discharge  mode  for  two 
types  of  MEAs  match  quite  satisfactorily,  showing  that  the 
electrodes  obtained  from  TVN  systems  have  consistent  cat¬ 
alytic  properties.  The  overall  performances  are  identical  in 
Fig.  5  for  both  the  segmented  and  unsegmented  MEAs,  irre¬ 
spective  of  the  discharge  mode.  Based  on  the  results  seen  in 
Fig.  5  it  can  be  inferred  that  there  are  no  significant  differ¬ 
ences  in  terms  of  overall  activation,  ohmic  and  mass  trans¬ 
port  properties  between  the  two  types  of  MEA  and  that  the 
discharge  mode  does  not  affect  the  overall  cell  performance. 
Slightly  greater  area  between  the  current  collectors  in  the 
common  MEA  seems  to  have  a  negligible  effect  probably 
due  to  lack  of  electronic  contact. 

Having  established  that  the  overall  properties  of  the 
two  MEAs  are  quite  similar,  the  question  is  whether  it  is 
significant  local  variations  in  the  ohmic,  transport  or  ki¬ 
netic  properties  along  the  length  of  the  channel  that  causes 
non-uniformity  in  the  case  of  common  MEA  under  po¬ 
tentiostatic  discharge.  Ideally,  if  all  relevant  properties  are 
invariant  along  the  channel  length,  the  individual  segment 
performances  should  be  identical  under  the  chosen  op¬ 
erating  conditions  (that  ensures  no  change  in  reactant  or 
product  concentrations  along  the  channel)  irrespective  of 


the  electrode  configuration  or  discharge  mode.  If  this  was 
possible,  segmenting  the  electrode  would  be  unnecessary. 
However,  the  individual  segment  polarization  curves  shown 
here  suggest  that  there  are  indeed  variations  in  segment 
performances  in  the  case  of  the  common  MEA  which  are 
mitigated  by  segmenting  the  electrodes.  For  the  common 
MEA,  the  fact  that  discharge  mode  has  an  effect  on  the 
variations  in  segment  performances  suggest  some  kind  of 
segment  to  segment  interactions. 

In  a  segmented  MEA,  by  isolating  the  electrodes  physi¬ 
cally,  all  interactions  between  segments  are  eliminated  bar¬ 
ring  the  common  membrane  and  a  shared  gas  channel.  Hence 
one  would  expect  to  see  the  segments  act  as  independent 
single  cells  with  no  interference  from  the  neighboring  seg¬ 
ments  under  certain  operating  conditions  that  do  not  induce 
any  variations  in  reactant  and  product  concentrations.  How¬ 
ever,  variability  could  exist  between  segments  due  to  inher¬ 
ent  differences  in  local  membrane  properties  or  variations  in 
contact  resistances.  Under  subtle  operating  conditions  used 
here,  such  as  using  pure  humidified  hydrogen  and  oxygen 
and  low  reaction  demand  that  do  not  induce  significant  con¬ 
centration  variations  along  the  channel,  the  results  from  the 
segmented  MEA  suggest  no  interaction  between  segments 
as  expected  and  it  is  also  encouraging  to  note  that  segmented 
MEAs  can  be  fabricated  with  comparable  segment  proper¬ 
ties. 

In  a  common  electrode  there  can  be  interferences  be¬ 
tween  segments  due  to  differences  in  local  properties  of  the 
electrode.  Liquid  water  dynamics  between  the  regions  un¬ 
der  two  adjacent  current  collectors  may  play  a  role  in  these 
segment  to  segment  interactions.  However,  the  extreme  as¬ 
pect  ratio  between  the  direction  along  the  channel  and  that 
normal  to  the  reactive  interface  and  the  drastic  differences  in 
the  available  cross-sectional  area  for  flow  in  these  two  direc¬ 
tions  are  expected  to  significantly  limit  this  liquid  transport 
interaction.  Also,  liquid  water  flooding  effects  also  affect 
the  overall  performance  of  the  MEA  which  was  not  the  case 
here.  Moreover,  the  fact  that  the  common  MEA  shows  very 
similar  segment  performances  under  galvanostatic  mode 
and  not  under  potentiostatic  mode  (while  segmented  MEA 
results  were  similar  for  both  discharge  modes)  lead  us  to 
believe  that  the  interactions  could  be  electronic  and  not 
kinetic  in  nature.  Electronic  interactions  arise  mainly  due 
to  differences  in  contact  resistances  between  the  various 
current  collector  segments  and  the  GDL  (assuming  unifor¬ 
mity  in  GDL/electrode  electronic  properties).  However,  the 
similarity  between  the  two  types  of  discharges  and  lack  of 
significant  segment  to  segment  variations  in  the  segmented 
MEA  suggests  reasonable  uniformity  in  contact,  which  is 
counter-intuitive  to  our  previous  hypothesis  on  electronic 
interactions.  In-depth  information  on  the  potential  and  cur¬ 
rent  vector  distribution  within  the  diffusion  layer  for  the  two 
types  of  MEAs  under  different  discharge  modes  is  essential 
to  verify  the  validity  of  our  hypothesis  about  electronic 
interactions  in  a  common  MEA  and  to  obtain  a  clear  under¬ 
standing  of  the  phenomena  involved.  Hence  a  simple  model 


D.  Natarajan,  T.  Van  Nguyen/ Journal  of  Power  Sources  135  (2004)  95-109 


101 


Fig.  6.  Schematic  of  the  model  domain  that  includes  the  GDL,  one  half  of  two  adjacent  electrodes  separated  by  plastic  shoulder. 


was  developed  to  assist  in  the  interpretation  of  the  experi¬ 
mental  results  and  is  presented  in  the  following  sections. 

4.1.  Modeling 


At  the  boundary  3-B  : 


GDL 


( 


A0b 


\  ^CON,B 


A  two-dimensional  mathematical  model  was  developed 
to  map  the  solid  phase  potential  (0)  distribution.  The  model 
domain  is  presented  in  Fig.  6  along  with  the  various  bound¬ 
aries  and  dimensions.  The  domain  included  half  of  two  ad¬ 
jacent  current  collectors  in  contact  with  a  common  diffusion 
layer,  with  a  plastic  shoulder  separating  the  current  collec¬ 
tors.  The  model  accounts  for  the  dimension  in  the  chan¬ 
nel  direction  along  which  the  current  collector  segments  are 
placed  (X-direction)  and  the  dimension  normal  to  the  re¬ 
active  interface  (T-direction).  The  governing  equations  and 
boundary  conditions  are  provided  below. 

•  Governing  equations 

Within  the  gas  diffusion  layer  : 

/  d20  d20\ 

~  I  ax, GDL  7^2  +  ctKGDL^2  )  =  0  ^ 


Within  current  collectors  A  and  B 
fd20  d2<P  . 

-accl9^+9^1=° 


(2) 


Within  the  plastic  separator  :  0  =  0.0 


(3) 


•  Boundary  conditions 

At  boundaries  2,  4,  3-S,  5-A  and  5-B  :  n  •  V<£>  =  0 
where  n  is  the  unit  vector  in  the  normal  direction. 


At  the  boundary  3- A  : 


y&Y,  GDL 


/  IA^aI 

\^CON,A 


At  the  boundary  6- A  and  6-B  : 
30 


&cc 


dY 


=  T\pp  or  0  =  Lapp 


At  the  reactive  interface  1  : 

/  30\ 

-^KGDL-j 

=  -i0(lO[(d>_£°)/’,c]  -  io[-(d>-£°)/’'c])  (7) 

Here  crx,  gdl  and  crY, gdl  stand  for  the  in-plane, 
through-plane  conductivities  of  the  gas  diffusion  layer, 
c>cc  represents  the  bulk  conductivity  of  the  current  col¬ 
lector,  Rqon,a  and  7?con,b  are  the  areal  contact  resis¬ 
tances  between  current  collectors  A  and  B  and  the  GDL, 
and  A0a  and  A 0%  are  the  potential  drops  at  the  corre¬ 
sponding  interfaces  due  to  contact  resistance.  The  reader 
is  referred  to  the  ‘list  of  symbols’  section  for  explana¬ 
tion  on  other  variables.  The  partial  differential  equations 
were  discretized  using  finite  difference  method  and  the 
resulting  algebraic  equations  were  solved  using  a  banded 
matrix  solver  [28]. 

As  can  be  seen  from  the  boundary  conditions,  the  model 
was  setup  to  accommodate  both  galvanostatic  and  poten- 
tiostatic  simulations.  The  boundary  condition  at  interface 
1  was  set  to  represent  the  slow  oxygen  reduction  reaction 
in  most  cases.  The  model  was  used  to  simulate  galvano¬ 
static  and  potentiostatic  discharges  employing  segmented 
and  unsegmented  SIGRACET®  and  Toray®  gas  diffusion 
layers.  Toray®  GDL  with  different  electronic  properties 
was  included  in  these  simulations  to  validate  the  qualitative 
predictions  of  this  model.  The  relevant  model  parameters 
are  provided  in  Table  1.  The  exchange  current  density  was 
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Table  1 

Parameter  values  used  in  the  model 


SIGRACET®  GDL 

Toray®  GDL 

In-plane  conductivity  for  common  MEA,  ax,  gdl  (S/cm) 

78.5 

200.0 

In-plane  conductivity  for  segmented  MEA,  ax, gdl  (S/cm) 

0.785 

2.0 

Through-plane  conductivity,  a y  ,gdl  (S/cm) 

3.14 

14.0 

Areal  contact  resistance  of  collector  A,  Rcon.a  (^  cm2) 

0.0105 

0.0105 

Areal  contact  resistance  of  collector  B,  Rcon.b  (£2  cm2) 

0.021 

0.021 

Bulk  conductivity  of  collectors  A  and  B,  acc  (S/cm) 

670 

670 

Exchange  current  density  [29,30],  z'o  at  30  °C  (A/cm2) 

5.0  x  10-9 

2.5  x  10"9 

Equilibrium  potential,  (V) 

1.2 

1.2 

Cathode  Tafel  constant  [29,30],  ij c  at  30 °C  (V/decade) 

0.0693 

0.0693 

a  Parthasarathy  et  al.  [29]  reported  a  rj c  value  of  60  mV/decade  at  25  °C,  while  a  value  of  65  mV/decade  at  30  °C  was  suggested  by  Paik  et  al.  [30]. 


adjusted  to  provide  the  current  densities  in  the  range  of  the 
experimental  values.  The  exchange  current  density  values 
used  in  these  simulations  are  comparable  to  values  reported 
in  literature  for  ORR  kinetics  at  Platinum/Nafion®  interface 
at  30  °C  [29,30].  The  Tafel  slope  for  ORR  kinetics  chosen  for 
these  simulations  was  also  close  to  reported  literature  values 
[29,30].  The  through-plane  conductivity  of  SIGRACET® 
GDL  was  obtained  from  the  manufacturer.  The  in-plane 
conductivity  was  estimated  by  a  four  probe  conductivity 
measurement  technique  and  the  order  of  magnitude  of  the 
measured  value  was  also  verified  by  the  manufacturer.  Both 
in-plane  and  through-plane  conductivity  for  Toray®  gas  dif¬ 
fusion  layers  were  obtained  from  the  manufacturer  and  ver¬ 
ified  using  a  four  probe  milliohmmeter.  The  bulk  conductiv¬ 
ity  of  POCO  graphite  was  obtained  from  the  literature  [31]. 

The  contact  resistances  between  the  segmented  current 
collectors  and  the  diffusion  layers  were  estimated  experi¬ 
mentally.  The  segmented  plates  were  assembled  with  just 
a  SIGRACET®  or  Toray®  GDL  of  known  dimensions  un¬ 
der  a  compression  pressure  of  65  psi.  Moreover,  a  silicone 
gasket  was  used  as  a  frame  to  position  the  carbon  paper  ac¬ 
curately  on  the  segmented  current  collectors.  Though  these 
measurements  may  not  reflect  the  true  contact  resistances 
when  an  actual  MEA  is  used,  they  should  still  provide  a  rea¬ 
sonable  range.  This  type  of  assembly  includes  the  inherent 
through-plane  resistance  of  the  GDL  and  contact  resistances 
at  both  faces  of  the  GDL.  Measurements  were  conducted  on 
multiple  specimens  for  both  types  of  GDL  for  all  the  cur¬ 
rent  collectors  at  varying  current  densities  and  the  results 
are  tabulated  in  Table  2.  The  values  provided  in  Table  2 
were  based  on  the  total  geometric  area  of  the  GDL  and  not 
just  the  shoulder  area  of  the  current  collectors.  Hence  the 
actual  contact  resistances  encountered  is  probably  slightly 
lower.  The  range  of  the  measured  contact  resistance  values 


for  the  two  types  of  GDL  were  very  similar  as  expected.  To 
this  author’s  knowledge,  even  the  high  end  of  the  range  of 
contact  resistances  in  Table  2  is  among  the  lowest  reported 
in  the  literature.  Mench  and  Wang  [25]  reported  an  average 
contact  resistance  of  0.0417  £2  cm2.  Barbir  et  al.  [32]  sug¬ 
gest  contact  resistances  as  high  as  0.15  Q  cm2.  Mepsted  and 
Moore  [31]  conducted  a  comprehensive  study  on  different 
materials  used  for  bipolar  plates  and  concluded  that  graphite 
plates  offer  the  lowest  contact  resistance  (or  surface  resistiv¬ 
ity)  of  all  the  materials.  Their  measured  graphite-graphite 
contact  resistance  was  about  0.014  £2  cm2  (for  the  sake  of 
comparison).  Based  on  the  values  reported  in  Table  2  it  is 
also  seen  that  the  lower  end  of  the  range  of  areal  contact 
resistances  is  almost  the  same  as  that  of  the  through-plane 
areal  resistance  of  the  SIGRACET®  GDL. 

Initially  the  areal  contact  resistances  for  both  the  cur¬ 
rent  collectors  simulated  in  the  model  were  assigned  the 
same  value  (0.0105  £2  cm2)  as  that  of  the  through-plane  areal 
resistance  of  the  gas  diffusion  layer  (SIGRACET®).  The 
T-direction  thickness  of  the  current  collector  was  arbitrar¬ 
ily  chosen  to  be  330  [xm.  The  thickness  of  the  bulk  current 
collectors  is  not  crucial  as  the  potential  drop  in  the  current 
collectors  were  expected  to  be  negligible.  When  the  con¬ 
tact  resistances  were  uniform,  the  simulation  results  for  both 
discharge  modes  were  symmetrical  about  the  mid  point  on 
the  plastic  separator  and  the  current  or  voltage  responses 
at  the  two  collectors  were  identical  as  expected  (results  are 
not  provided  here  for  the  sake  of  brevity).  The  areal  con¬ 
tact  resistance  at  current  collector  B  was  then  doubled  to 
correspond  to  the  highest  measured  areal  contact  resistance 
(0.021  Q  cm2). 

The  results  of  this  simulation  (both  constant  current  and 
constant  voltage  mode)  for  a  SIGRACET®  GDL  are  pro¬ 
vided  in  Lig.  7.  Qualitatively,  the  model  predictions  match 


Table  2 

Contact  resistances  measured  for  SIGRACET®  GDL  and  Toray®  GDL 

Through-plane  areal  resistance  Range  of  measured  total  through-plane  Range  of  calculated  through-plane 

of  GDL  (£2  cm2)  (obtained  areal  resistance  (£2  cm2)  (includes  areal  contact  resistance  (£2  cm2) 

from  manufacturer)  GDL  and  contact  resistances) 

SIGRACET®  GDL  0.0105  0.032-0.046  0.011-0.018 

Toray®  GDL  0.0025  0.024-0.047  0.011-0.022 
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Abs.  value  of  segment  current  density  (A/cm2) 


Fig.  7.  Galvanostatic  and  potentio static  model  simulations  for  a  cathode 
made  of  common  SIGRACET®  GDL  with  non-uniform  contact  resistances 
(10.5  and  21m£2cm2  for  current  collectors  A  and  B,  respectively). 

quite  well  with  the  experimental  results  presented  in  Fig.  4. 
The  significant  discrepancy  in  the  ordinate  (cell  voltage) 
axis  can  be  attributed  to  the  fact  that  the  simple  model  de¬ 
veloped  here  does  not  account  for  the  potential  losses  at  the 
membrane  and  anode.  Nevertheless  the  model  simulations 
capture  the  differences  encountered  between  the  two  types 
of  discharge  modes.  It  is  evident  that  the  differences  in  the 
polarization  behavior  between  constant  current  and  constant 
voltage  mode  observed  in  experiments  is  clearly  due  to  dif¬ 
ferences  in  contact  resistances  between  the  current  collector 
segments  and  the  GDL.  It  should  be  noted  that  the  kinetic 
parameters  were  held  constant  in  these  simulations.  The  case 
of  segmented  electrodes  was  also  simulated  by  arbitrarily  de¬ 
creasing  the  in-plane  conductivity  of  the  GDL  by  two  orders 
of  magnitude.  The  results  are  presented  in  Fig.  8.  Similar  to 
the  experimental  data  shown  in  Fig.  3,  when  the  electrodes 
are  segmented,  the  polarization  curves  lie  on  top  of  each 
other  irrespective  of  the  discharge  mode.  One  of  the  major 
concerns  among  the  fuel  cell  research  community  regard¬ 
ing  segmenting  electrodes  is  achieving  uniform  comparable 
contact  and  performance  from  segment  to  segment  [21-27]. 


Fig.  8.  Galvanostatic  and  potentio  static  model  simulations  for  a  cathode 
with  segmented  SIGRACET®  GDL. 


The  experimental  and  model  results  suggest  that,  small  dif¬ 
ferences  in  contact  resistances  do  not  significantly  manifest 
themselves  when  the  electrode  is  segmented.  On  the  other 
hand  despite  similarity  in  kinetic  and  transport  properties 
the  same  difference  in  contact  resistances  significantly  af¬ 
fects  the  results  in  the  case  of  the  common  MEA,  especially 
in  the  constant  voltage  mode. 

To  further  understand  this  behavior,  the  potential  and  elec¬ 
tronic  flux  vector  distributions  within  the  common  GDL  are 
mapped  in  Figs.  9  and  10  for  the  two  types  of  discharges. 
Figs.  9a  and  10a  clearly  show  a  symmetric  potential  and 
electronic  flux  distribution  within  the  gas  diffusion  layer  for 
the  case  of  the  galvanostatic  simulation.  The  solid  potential 
has  a  greater  drop  at  the  interface  between  the  current  collec¬ 
tor  B  and  the  GDL  due  to  the  higher  contact  resistance.  As 
expected,  the  potential  drop  within  the  current  collectors  is 
negligible.  The  potentiostatic  simulations  shown  in  Figs.  9b 
and  10b  show  a  drop  in  the  potential  within  the  GDL  from 
the  region  above  current  collector  B  towards  the  region  over 
current  collector  A  indicating  the  existence  of  in-plane  cur¬ 
rent  flow  within  the  GDL.  The  electronic  flux  vector  plot 
confirms  this  phenomenon  indicating  flow  of  electrons  from 
current  collector  A  towards  the  regions  over  current  collec¬ 
tor  B.  Thus  despite  the  significant  geometric  size  differences 
associated  with  the  in-plane  and  through-plane  directions, 
there  exists  significant  interactions  between  segments.  This 
is  probably  due  to  the  inherent  higher  in-plane  conductivity 
compared  to  through-plane  values  of  typical  GDL  material 
(see  Table  1)  arising  from  carbon  fiber  orientation.  Hottinen 
et  al.  [27]  used  the  through-plane  resistance  of  the  GDL  in 
all  directions  in  their  model  and  hence  did  not  observe  a 
similar  phenomenon. 

Ignoring  bulk  collector  losses  the  total  voltage  drop 
between  the  reactive  interface  and  the  voltage  sensing 
point  in  the  current  collector  can  be  expressed  as  Af  = 
AVcon  +  AVgdl  +  AVk,  where  AVcon  is  the  potential 
loss  associated  with  contact,  AVgdl  is  the  loss  across 
the  diffusion  layer  and  AVk  is  the  activation  loss.  During 
galvanostatic  control,  AVgdl  and  AVk  are  fixed  by  the 
current  demand  per  segment.  Hence  the  only  difference 
between  the  segments  would  arise  out  of  the  differences 
in  contact  between  current  collectors  and  GDL,  which  are 
quite  small.  In  the  case  of  potentiostatic  control,  however, 
all  the  three  losses  mentioned  above  vary.  The  potential 
at  the  reactive  interface  depends  on  both  the  contact  and 
GDL  losses  and  due  to  the  exponential  nature  of  the  ki¬ 
netic  term,  differences  in  the  solid  potentials  at  the  reactive 
interface  can  result  in  significant  variations  in  the  current 
distribution. 

Fuel  cell  researchers  usually  assume  uniformity  in  local 
current  density  distribution  in  the  kinetic  region.  However, 
these  results  clearly  indicate  that  in  unsegmented  electrodes, 
non-uniformity  in  contact  resistances  manifest  themselves 
significantly  (depending  on  the  discharge  mode),  even  in 
the  low-current  density  activation  region  and  such  an  as¬ 
sumption  is  not  appropriate  if  uniform  contact  cannot  be 
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Fig.  9.  Model  results  of  potential  distribution  within  the  SIGRACET®  GDL  of  the  cathode  for  (a)  galvanostatic  discharge  and  (b)  potentiostatic  discharge. 


achieved.  Though  this  study  is  limited  to  the  low  current 
density  region,  the  effects  of  contact  variations  will  persist 
over  the  entire  operating  range  of  potentials  (or  current 
density)  and  can  be  expected  to  grow  in  significance  with 
increasing  current  densities.  However,  transport  pertaining 
to  liquid  water  dynamics  will  eventually  start  to  gain  sig¬ 
nificance  and  this  phenomenon  might  mask  the  electronic 


interactions  between  the  segments.  Hence,  it  is  essential 
to  first  study  the  effect  of  contact  resistance  variations  at 
low  current  density  region  where  liquid  water  influences 
are  minimal.  Once  these  contact  issues  are  addressed,  the 
interactions  between  segments  at  higher  current  densities 
can  be  easily  discerned  from  liquid  water  effects,  given  the 
linear  dependence  of  ohmic  effects  on  current  densities. 
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(a)  Electrode  width  (cm) 
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(b)  tlectrode  width  (cm) 

Fig.  10.  Model  results  of  current  vector  distribution  within  the  SIGRACET®  GDL  of  the  cathode  for  (a)  galvanostatic  discharge  and  (b)  potentio static 
discharge  (arrow  sizes  not  scaled  to  the  magnitude  of  the  vector). 
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It  should  be  noted  that  the  experimental  results  were  gen¬ 
erated  using  either  segmented  or  unsegmented  electrodes  on 
both  the  anode  and  cathode.  Hence  differences  in  contact 
resistances  on  the  anode  side  were  also  evaluated  using  the 
model.  Eq.  (7)  that  represents  the  boundary  condition  at  the 
reactive  interface  was  replaced  by  a  simple  linear  expression 
shown  in  Eq.  (8)  to  capture  the  facile  anode  kinetics. 

d&\  0-E° 

°F,GDL—  I  =  -  (8) 

3  Y  )  r)  A 

where  ?/.\  has  a  value  [33]  of  40mV/A/cm2.  All  other  pa- 
rameters  were  kept  the  same  as  that  of  the  cathode  simula¬ 
tions.  Once  again  the  model  simulations  were  symmetrical 
about  the  plastic  separator  midpoint  and  the  segment  per¬ 
formances  were  identical  when  the  contact  resistances  were 


kept  the  same  for  the  two  current  collectors.  The  simula¬ 
tion  results  for  current  collectors  with  different  contact  re¬ 
sistances  are  presented  in  Fig.  11a  and  b  for  a  common  and 
segmented  anode.  These  results  are  significantly  different 
from  the  cathode.  While  the  linear  current-voltage  response 
is  expected,  it  can  be  seen  that  the  two  segments  do  not  per¬ 
form  uniformly  in  both  MEAs  irrespective  of  the  discharge 
mode.  Like  the  cathode,  during  potentiostatic  runs  on  the  an¬ 
ode,  all  the  potential  losses  associated  with  current  flow  are 
variables  and  are  ultimately  related  to  the  differential  con¬ 
tact,  leading  to  performance  variations  between  segments. 
In  the  case  of  galvanostatic  mode,  once  again  the  difference 
between  segment  performances  is  mainly  controlled  by  the 
differences  in  contact.  Unlike  the  cathode  where  the  large 
kinetic  losses  masks  these  differences,  in  the  anode  the  con¬ 
tact  losses  are  comparable  to  the  activation  loss  and  manifest 
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2 

(a)  Segment  current  density  (A/cm  ) 


Seg.# 


Solid  line  -  Galvanostatic,  Dashed  line  -  Potentiostatic  v  "  1 

□  -  2 

Fig.  11.  Galvanostatic  and  potentiostatic  model  simulations  for  an  anode 
with  (a)  common  SIGRACET®  GDL  and  (b)  segmented  SIGRACET® 
GDL. 


Fig.  12.  Galvanostatic  and  potentiostatic  segment  polarization  curves  using 
a  segmented  MEA  with  Toray®  GDL.  H2  flow  rate  =  2.3  A/cm2,  O2 
flow  rate  =  2.3  A/cm2,  H2  humidifier  temperature  =  70  °C,  O2  humidifier 
temperature  =  25  °C  and  cell  temperature  =  30  °C. 


Abs.  value  of  segment  current  density  (A/cm2) 


themselves  even  in  the  constant  current  mode  of  operation 
leading  to  some  interaction  between  segments.  One  impor¬ 
tant  point  to  note  is  that  the  differences  in  the  simulated  per¬ 
formance  from  segment  to  segment  is  quite  negligible  at  the 
anode  (ca.  1  mV)  when  compared  to  the  cathode  (10-12  mV) 
at  a  given  current  density  (ca.  100  mA).  Hence  the  cathode 
contributes  mostly  to  the  differences  between  segment  per¬ 
formances  seen  in  the  experiments.  These  model  results  for 
the  anode  suggest  that  any  difference  in  contact  will  show 
up  irrespective  of  the  whether  the  electrode  is  segmented  or 
unsegmented.  However,  the  differences  caused  by  the  dif¬ 
ferential  contact  at  the  anode  have  negligible  contribution  to 
the  overall  interaction  between  segments. 

Experiments  and  model  simulations  were  also  conducted 
on  Toray®  gas  diffusion  layers  to  further  verify  the  observa¬ 
tions  discussed  in  the  previous  sections.  The  Toray®  GDL 
has  a  higher  bulk  density  and  lower  porosity.  Hence  the  con¬ 
ductivities  in  both  the  in-plane  and  through-plane  directions 
are  much  higher  (see  Tables  1  and  2).  Figs.  12  and  13  pro¬ 
vide  the  experimental  data  on  segmented  and  unsegmented 
(common)  MEAs  prepared  using  Toray®  carbon  paper  while 
the  overall  performances  are  compared  in  Fig.  14.  The  trends 


Fig.  13.  Galvanostatic  and  potentiostatic  segment  polarization  curves  us¬ 
ing  a  common  MEA  with  Toray®  GDL.  H2  flow  rate  =  2  A/cm2,  O2 
flow  rate  =  2  A/cm2,  H2  humidifier  temperature  =  70  °C,  O2  humidifier 
temperature  =  25  °C  and  cell  temperature  =  30  °C. 


Fig.  14.  Comparison  of  total  cell  current  of  common  and  segmented 
MEAs  under  galvanostatic  and  potentiostatic  discharge  with  Toray®  GDL. 
H2  flow  rate  =  2  or  2.3  A/cm2,  O2  flow  rate  =  2  or  2.3  A/cm2,  H2 
humidifier  temperature  =  70  °C,  O2  humidifier  temperature  =  25  °C  and 
cell  temperature  =  30  °C. 
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Fig.  15.  Galvanostatic  and  potentiostatic  model  simulations  for  the  cathode 
with  a  common  Toray®  GDL. 

represented  in  these  figures  are  quite  similar  to  results  from 
SIGRACET®  GDL.  Here  again,  segmented  MEAs  provide 
polarization  curves  that  are  almost  identical  irrespective  of 
the  discharge  mode.  The  common  MEA  shows  significant 
differences  in  the  individual  collector  performances  under 
potentiostatic  mode  while  the  galvanostatic  curves  are  very 
close  to  each  other.  As  before,  Fig.  14  does  not  indicate 
any  significant  difference  between  the  two  MEAs  in  terms 
of  overall  performance  validating  the  need  for  a  segmented 
electrode  and  current  collector  setup  for  local  current  den¬ 
sity  distribution  studies.  The  individual  segment  and  overall 
performances  were  in  general  lower  than  the  SIGRACET® 
GDL.  This  was  probably  due  to  the  lower  porosity  of  the 
Toray®  diffusion  layer  resulting  in  poorer  gas  and  liquid  wa¬ 
ter  transport  leading  to  lower  oxygen  concentration  at  the 
interface.  Interestingly,  in  the  constant  voltage  mode,  the 
segments  (3  and  4)  that  tends  to  deviate  from  the  galvanos¬ 
tatic  curves  seems  to  be  consistent  between  the  two  type  of 
diffusion  layers  (SIGRACET®  and  Toray®). 

The  corresponding  model  results  for  common  and  seg¬ 
mented  MEAs  fabricated  with  Toray®  GDL  are  presented  in 
Figs.  15  and  16.  In  these  simulations,  the  effective  exchange 


Fig.  16.  Galvanostatic  and  potentiostatic  model  simulations  for  the  cathode 
with  a  segmented  Toray®  GDL. 


current  density  was  reduced  by  half  to  allow  for  the  model 
simulations  to  roughly  match  the  current  density  range  ob¬ 
tained  in  the  experiments.  Other  than  exchange  current  den¬ 
sity  and  the  relevant  conductivity  values,  no  other  changes 
were  made  to  the  model  parameters.  The  contact  resistances 
were  kept  the  same  as  before  since  experimentally  measured 
values  for  the  two  types  of  GDL  were  in  the  same  range. 
The  model  results  in  Figs.  15  and  16  qualitatively  match 
the  experimental  results  for  Toray®  GDL.  Once  again  the 
spread  in  the  polarization  curves  for  a  common  MEA  on 
constant  voltage  mode  caused  by  differential  contact  is  con¬ 
sistent  with  experimental  results. 

The  model  and  experimental  results  from  the  two  types 
of  diffusion  layers  clearly  show  that  inter- segment  interac¬ 
tions  in  common  MEA  can  be  avoided  only  if  the  differences 
in  segmented  current  collector- MEA  contact  resistances  are 
minimal.  As  mentioned  before,  the  contact  resistances  mea¬ 
sured  in  this  study  are  among  the  lowest  reported  values  and 
achieving  further  uniformity  at  this  scale  can  be  very  diffi¬ 
cult.  Moreover,  the  same  differences  in  contact  that  also  ex¬ 
isted  in  the  case  of  a  segmented  MEA  did  not  significantly 
affect  the  uniformity  in  the  segment  performances.  It  should 
also  be  noted  that  when  using  segmented  electrodes,  the  dis¬ 
tance  between  the  segments  (and  hence  the  segmented  cur¬ 
rent  collectors)  is  not  crucial  as  long  as  no  concentration 
variations  are  introduced  along  the  channel. 

Finally,  the  question  of  whether  the  common  membrane 
affects  the  current  density  distribution  in  terms  of  segment 
to  segment  interaction  needs  to  be  addressed.  The  ionic 
conductivity  of  even  a  well  hydrated  membrane  is  usually 
about  two  orders  of  magnitude  smaller  than  the  electronic 
conductivity  of  the  electrode.  Though  the  local  membrane 
conductivity  is  a  strong  function  of  the  hydration  state  of 
the  membrane  which  can  vary  along  the  channel  depend¬ 
ing  on  the  operating  conditions  like  flow  rates,  humidity, 
temperature,  etc.,  currently  there  is  no  evidence  of  any 
anisotropy  in  conductivity  of  the  membrane  used  in  this 
study,  unlike  the  electrode.  In  other  words,  the  overall  con¬ 
ductivity  of  the  membrane  is  expected  to  be  the  same  in  the 
normal  direction  and  the  direction  along  the  channel.  Also, 
typical  aspect  ratios  related  to  the  membranes  (50  jxm  by 
few  mm  for  Nafion®  112)  are  also  more  extreme  than  that 
of  the  electrode.  (175-350  jam  by  few  mm).  Based  on  these 
observations,  minimal  segment  to  segment  interactions  or 
‘crosstalk’  is  expected  across  the  membrane  when  a  MEA 
with  segmented  electrodes  and  common  membrane  is  used 
for  current  density  distribution  studies.  Qualitative  evidence 
for  this  hypothesis  was  presented  by  Natarajan  and  Nguyen 
[34].  The  authors  showed  that  the  differences  in  segment 
performances  were  insignificant  for  segmented  electrode 
MEAs  with  common  and  segmented  membrane. 

As  mentioned  earlier  in  the  introduction  section,  sophisti¬ 
cated  models  have  been  developed  to  capture  the  two-phase 
transport  and  electrokinetic  phenomena  in  a  PEM  fuel  cell. 
Reasonably  accurate  estimates  of  the  various  physical  and 
chemical  properties  of  the  components  of  the  PEM  fuel  cell 
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are  essential  to  transform  these  models  into  predictive  tools, 
for  which  relevant  experimental  data  are  needed.  Important 
kinetic  parameters  such  as  the  intrinsic  exchange  current 
densities  for  platinum  catalyzed  H2  oxidation  and  O2  re¬ 
duction  reactions  have  been  experimentally  estimated  and 
methods  to  extrapolate  these  results  to  porous  catalyst  layers 
such  as  in  PEM  fuel  cells  have  been  established  [29,30,35]. 
However,  there  is  a  lack  of  information  on  physical  proper¬ 
ties  of  the  gas  diffusion  layers  such  as  gas  and  liquid  water 
permeability  and  their  dependence  on  liquid  water  satura¬ 
tion,  the  functional  dependence  of  capillary  pressure  on  sat¬ 
uration,  etc.  Fuel  cell  researchers  have  attempted  to  come 
up  with  estimates  for  these  gas  and  liquid  transport  param¬ 
eters  by  comparing  experimental  data  to  model  predictions 
[14,17,36].  Natarajan  and  Nguyen  [14]  for  instance,  stud¬ 
ied  the  effect  of  model  domain  selection  (2D  versus  3D)  on 
permeability  and  capillary  function  estimates  for  the  GDL 
in  a  PEM  fuel  cell.  Their  work,  suggests  a  clear  need  for 
a  robust  three  dimensional  model  along  with  experimental 
data  on  local  current  density  distribution  at  various  operat¬ 
ing  conditions  to  accurately  estimate  such  GDL  properties. 
Segmented  fuel  cell  assemblies  are  needed  to  generate  such 
local  current  density  distribution  data. 

The  information  provided  in  this  paper  clearly  demon¬ 
strates  that  the  significant  difference  between  the  in-plane 
and  through-plane  conductivities  of  the  GDL  and  variations 
in  contact  resistances  affects  the  local  current  density  distri¬ 
bution,  when  only  the  current  collector  is  segmented  and  the 
electrode  is  not.  Experimental  results  from  such  a  setup  can 
only  be  compared  to  models  where  the  electronic  properties 
of  the  GDL  and  contact  resistance  variations  are  properly 
accounted  for.  Most  of  the  models  that  is  available  in  the 
literature  deal  with  an  unsegmented  electrode  and  do  not  ex¬ 
plicitly  account  for  the  electronic  properties  of  the  GDL  and 
variations  in  contact  resistance  which  affects  the  solid  poten¬ 
tial  distribution  and  hence  the  local  current  density.  A  few 
models  [9]  do  account  for  the  electronic  conduction  in  the 
solid  phase,  but  only  in  terms  of  a  heat  source  (Joule  heating) 
that  affects  the  temperature  distribution.  Such  an  approach 
can  only  be  justified  if  it  is  assumed  that  the  effect  of  varia¬ 
tions  in  contact  resistance  is  minimal.  The  results  provided 
here  clearly  show  that  by  segmenting  the  electrode  this  as¬ 
sumption  of  minimal  contact  effects  is  satisfied.  Then  the 
remaining  questions  are:  (1)  whether  there  exists  significant 
liquid  water  movement  within  the  gas  diffusion  layer  in  the 
direction  along  the  channel  which  might  affect  the  local  cur¬ 
rent  density  distribution  (i.e.  along  the  direction  suggested 
for  electrode  segmentation)  and  would  segmenting  the  elec¬ 
trode  unrealistically  eliminate  such  a  transport  feature?  and 
consequently  (2)  can  data  from  segmented  electrode  studies 
be  compared  to  model  predictions  based  on  common  elec¬ 
trodes?  As  mentioned  before,  owing  to  the  extreme  aspect  ra¬ 
tio  of  the  PEM  fuel  cell  electrode  and  the  related  differences 
in  area  for  flow,  it  is  expected  that  liquid  water  transport 
within  the  GDL  along  the  channel  will  be  minimal  and  seg¬ 
menting  the  electrodes  would  still  capture  the  current  density 


distribution  that  one  can  expect  in  an  ideal  (with  minimal 
contact  influences)  unsegmented  electrode.  This  hypothesis 
can  be  further  justified  by  the  fact  that  while  transport  prop¬ 
erties  within  the  GDL  are  expected  to  be  isotropic  (as  long 
as  the  structure,  morphology  and  Teflon  content  distribu¬ 
tions  are  uniform  across  the  GDL),  current  GDL  materials 
do  exhibit  significant  anisotropy  when  it  comes  to  electronic 
conductivity. 

5.  Conclusions 

Current  density  distribution  measurements  under  gal- 
vanostatic  (constant  current)  and  potentio static  (constant 
voltage)  modes  were  carried  out  on  membrane  electrode 
assemblies  fabricated  with  segmented  and  unsegmented 
electrodes  using  a  segmented  current  collector  and  flow  field 
setup.  Two  commercially  available  gas  diffusion  layers, 
namely  SIGRACET®  GDL  30  BC  from  SGL-CARBON 
Inc.  and  Toray®  carbon  paper  TGPH-120  from  ETEK  Inc., 
were  studied. 

For  both  types  of  GDL,  when  the  electrode  was  not  seg¬ 
mented,  the  individual  segment  performances  showed  sig¬ 
nificant  variations  when  discharged  at  constant  voltage  and 
the  variations  were  minimal  under  galvanostatic  mode.  The 
overall  cell  performances,  however,  were  very  similar  be¬ 
tween  the  two  types  of  discharges. 

In  the  case  of  both  the  SIGRACET®  and  Toray®,  the  po¬ 
larization  curves  of  the  individual  segments  did  not  show 
significant  differences  irrespective  of  the  discharge  mode 
when  the  electrodes  were  segmented.  Here  again  the  overall 
performances  of  the  MEA  were  identical  between  the  con¬ 
stant  current  and  constant  voltage  modes. 

A  simple  two-dimensional  mathematical  model  was  de¬ 
veloped  to  map  the  potential  distribution  in  the  electrode. 
The  model  simulations  qualitatively  matched  the  experimen¬ 
tal  results  and  also  proved  that  variation  in  contact  resis¬ 
tances  between  the  current  collectors  and  GDL  especially 
on  the  cathode  side  was  the  primary  cause  for  phenomenon 
observed  in  the  experiments. 

Based  on  the  experimental  and  model  results,  for  current 
density  distribution  studies,  segmenting  the  electrode  along 
with  the  current  collectors  is  recommended.  Inherent  small 
variations  in  contact  resistances  between  the  various  current 
collectors  and  GDL  do  not  significantly  manifest  them¬ 
selves  in  the  case  of  segmented  electrode  configuration, 
irrespective  of  discharge  mode.  On  the  contrary,  when  using 
a  common  electrode,  the  same  inherent  contact  variations 
will  significantly  affect  the  current  density  distribution  when 
the  segments  are  discharged  in  the  potentiostatic  mode.  It  is 
preferable  to  employ  the  constant  current  mode  of  discharge 
for  common  electrode  configurations  to  minimize  variations 
in  the  performance  due  to  contact  differences.  However,  it 
should  be  recognized  that  applying  constant  current  density 
at  every  current  collector  does  not  represent  actual  fuel 
cell  operating  condition  where  only  the  average  current 
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density  over  the  entire  fuel  cell  (single  cell  or  stack)  is 
specified. 
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